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Abstract 

We construct the evolution equations for the twist-3 chiral-odd spin independent fracture 
functions in QCD. The Gribov-Lipatov reciprocity relation is fulfilled at the one-loop level for 
the quasi-partonic two-particle cut vertices only. It is found that the rang of the anomalous 
dimensions matrix is infinite for any given moment of the three-parton fracture function as 
distinguished from the case of DIS-distributions where the rang of the matrix was finite and 
increases with the number of the moment. In the multicolour limit N c — > oo the evolution 
equation for the quark-gluon-quark correlation function decouples from another equation in 
the system and becomes homogeneous provided we discard the quark mass effects. This fact 
provides an opportunity to find its analytic solution explicitly in nonlocal form similarly to the 
DIS. 



1 Introduction. 



The deep inelastic scattering (DIS) of leptons on the hadron target is the most effective ex- 
perimental tool for studying the dynamics of hadron reaction on the parton level which has 
a firm basis in the quantum field theory provided by light-cone Operator Product Expansion 
(OPE) [J]]. The latter gives a strict theoretical ground for separation of two different scales of 
the underlying process: nonperturbative information is concentrated in the tower of local com- 
posite operators while the coefficient function which characterizes the hard interaction process 
of constituents can be dealt perturbatively. It makes possible the study of logarithmic violation 
of Bjorken scaling as well as the power suppressed contributions (higher twists) responsible for 
many subtle phenomena in a polarized scattering. However, there exists an equivalent approach 
to the analysis of the corresponding quantities which is based on factorization theorems |J and 
the evolution equations || J|]. In spite of the fact that the latter approach has some difficulties 
as compared to the former in the study of higher twists, like the loose of the explicit gauge and 
Lorentz invariance of calculations and also the presence of the overcomplete set of correlation 
functions, it has the important advantage as being the most close to the intuitive physical pic- 
ture and the similarity to the parton model. There is another advantage of the latter approach 
for studying of the higher twist effects from the point of view of experimental capabilities since 
the OPE provides us for the moments of the structure functions and in order to extract the 
former one needs to measure the latter in the whole region of the momentum fraction very 
accurately. Obviously, it is quite difficult task even for next generation colliders. While, with 
a set of evolution equations at hand, one can find, in principle, the Q 2 -dependence of the cross 
section in question by putting the experimental cuts on the region of the attained momentum 
fractions. This approach can also be used in the situations when the OPE is no longer valid. 
These are the inclusive production of the hadron in the e + e~-annihilation, semi-inclusive deep 
inelastic scattering, Drell-Yan lepton production et ctr. 

There is continuous interest in the inclusive production of hadron in hard reaction. This 
process involves a quark fragmentation function to describe the hadron production from the 
underlying hard parton scattering. But it differs considerably from the DIS as the short distance 
expansion is no longer relevant although the given process goes near the light cone. The 
theoretical basis for strict analyses of the above phenomena is provided by the generalization of 
the OPE to the time- like region in terms of Mueller's £-space cut vertices |5], which moments 
are essentially nonlocal in the coordinate space. Again this approach has all attractive features 
of the OPE as it provides a consistent framework to account for the higher twist effects 
as well as it allows to sum up the UV logarithms || by using the powerful methods of the 
renormalization group. 
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The semi-inclusive hadron production from a quark fragmentation is described in QCD by 
the specific nonperturbative correlation functions of quark and gluon field operators over the 
hadron states which can be identified with £-space cut vertices. While the behaviour of the 
latter with respect to the fraction of the parton momentum carried by the hadron is determined 
by the nonperturbative strong interaction dynamics, the large Q 2 -scale dependence is governed 
by perturbation theory only. 

In this paper we begin our study of the evolution of the chiral-odd and chiral-even twist-3 
fragmentation functions. We begin with nonpolarized twist-3 fracture functions. This is the 
simplest case with respect to the number of correlation functions involved in mixing under 
renormalization group evolution 1 . From the phenomenological point of view they appear, for 
example, in the cross section for semi-inclusive hadron production in the process of measuring 



the nucleon's transversity distribution h\(x) through deep inelastic scattering fLl 
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here « = 1 + Ax 2 M 2 /Q 2 , y = 1 - E'/E and the cross section is expressed in a frame where 
the lepton beam with energy E defines the z-axis and the x — z-plane contains the nucleon 
polarization vector, which has the polar angle x an d the scattered electron E' has the polar 
angles 9, <p. The functions G± and Gt are expressed in terms of the product of the distribution 
and fragmentation functions in the following way 

Gi(x,o = Iy:qm(x)f(o, (2) 



^(x)^(C) + ^l(x)r(C) 



(3) 

All of them have the expressions in QCD in terms of the light-cone Fourier transformation of 



correlation functions of fundamental quark and gluon fields over specific hadron states |TT|, [12 

I rdX 
1 rdX 



S+9i(x) = - / — e iA *W(0) 7+ 75^(An)|AO, 
I J Ztt 



SMx) = - I -e^(N\m^ + ^M^)\N), 



x For completeness we should note note that the Q 2 -evolution of the chiral-odd nonpolarized and polarized 
distributions was investigated in ref. || in the framework of the OPE by studying the mixing of composite 
local operators under renormalization, while the coordinate space evolution equations in multicolour QCD were 
constructed and solved in ref. Jl(| 
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HO = ±J^(0\ 1+ i>(\n)\h,X)(h,X\m\0)- (4) 

and X is given by equation (0). The summing over final state hadrons in the definition of the 
fracture functions is implicit throughout the paper. Note, that the physical regions are different 
for distribution and fragmentation functions: < x < 1 and 1 < ( < oo, respectively. 

In the present study we address ourselves to the problem of construction the evolution 
equations for the function I which mixes with other cut vertices during the renormalization. 
The outline of the paper is the following. In section 2 we give the definitions of the £-space cut 
vertices which are closed with respect to the renormalization group evolution and discuss the 
advantages of the light-cone gauge. In the third section we show our technique on an example 
of abelian gauge theory: we review the renormalization of the theory in the axial gauge and 
construct corresponding equations. In the section 4 the Q 2 -evolution equations are derived for 
the QCD case. The final section is devoted to the discussion of our results and conclusions. In 
the appendix we present some useful formulae which make the discussion more transparent. 

2 Definitions. 

It is well known that in order to endow with parton-like interpretation of field theoretical 
quantities and to get much deeper insight into the corresponding perturbative calculations it is 
necessary to use to ghost-free gauges. Owing to this fact we chose in what follows the light-cone 
gauge B + = n^B^ = for the boson field and make the Sudakov decomposition of the parton 
four-momentum into transverse and longitudinal components 

fc" = zjf + an» + fc£, (5) 

where p and n are null vectors which pick out two different directions on the light cone such 
that p 2 = n 2 = 0, (pn) = 1. The advantage of this gauge is that the gluon field operator B p is 
related to the field strength tensor G pa via simple relation 

1 r°° 

B p (Xn) = d^G+^Xn) = - / dze(X - z)G +fl (z) J (6) 

so, that the gauge invariant result can be restored after all required calculations have been 
performed. 

For our purposes it is much more suitable to deal with correlation functions listed below 2 

J(C) = \S ^ lK (0mXn)\h,X)(h,X\mm, (7) 
2 See discussion of this point in refs. Q 
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M(0 = ^e^(0|m 7+ ^(An)|/ i ,X)(/ i ,X|^(0)|0), 
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lX ^\0\ gl ^(\n)\h,X)(h,X\mBUfin)\0) : 
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^(CO) 
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4C7 ^ge^-^'(0|^(An)|/ i ,X)(/.,X|^(/in) g7 ^(0)|0). 



(9) 
(10) 

(11) 
(12) 



The quantities determined by these equations form the closed set under the renormalization, 
however, they are not independent since there is relation between them due to equation of 
motion for the Heisenberg fermion field operator 



i(0-M(0 + 1 <^i(C',C) = o. 

Here and in the following discussion we introduce the convention 



(13) 



(14) 



While the former two functions 1 and M. can be made explicitly gauge invariant by inserting 
the P-ordered exponential (which is unity in the gauge we have chosen) between the quark 
fields, the latter can be written in the gauge invariant way introducing the following objects: 



7^ (CO = C'-ZiCCC), ^2 (CO = CZ 2 (C,C') 



(15) 



Taking into account eq. (^) it is easy to verify that they are indeed expressed in terms of 
correlators involving the gluon field strength tensor. The functions ZS^ and Z^ are related 
by complex conjugation 



r{2) 



(1), 



'(2), 



(16) 



Their support properties can be found by applying Jaffe's recipe ||15|| . It has been shown that 
the field operators entering the definition of the correlation functions can be placed in the 
arbitrary order on the light cone with appropriate sign change according to their statistics. 
Then taking the particular ordering and saturating the correlation function by the complete 
set of the physical states we immediately obtain (for definiteness, we consider the function Z^) 

^i (1) (C',C) = i £ *(C - 1 - Cx)S(C - c Y )(o\mx)(h,x\i,\Y)(Y\B ± \o) 

4( > X,Y 

= h £ <KC - 1 - Cx)S(( -C- C Y )(0\^\h,X)(h,X\B ± \Y}(Y\m, (17) 
4( > XY 
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with (x, Cy > and we omit the unessential Dirac matrix structure of the vertex. From these 
equations the restrictions emerge on the allowed values of the momentum fractions: 1 < ( < oo, 
< £' < (. By analogy one can easily derive similar support properties for other functions. 



3 Construction of the evolution equations. 

Due to ultraviolet divergences of the momentum integrals in the perturbation theory there is 
logarithmic dependence of the parton densities on the normalization point. This dependence 
is governed by the renormalization group. The evolution equations for the leading twist cor- 
relation functions determining their Q 2 dependence can be interpreted in terms of the kinetic 
equilibrium of partons inside a hadron (for distribution functions) or hadrons inside a parton 
(for fragmentation function) under the variation of the ultraviolet transverse momentum cut-off 
0. However, beyond the leading twist the probabilistic picture is lost due to quantum me- 
chanical interference and more general quantities emerge, i.e. multiparticle parton correlation 
functions, whose scale dependence are determined by Faddeev type evolution equation with 



pairwise particle interaction |13|, 14 



There are two sources of the logarithmic dependence of the correlation functions. The 
first is the divergences of the transverse momentum integrals of the particles interacting with 
the vertex and forming the perturbative loop. Another source is the divergences due to the 
virtual radiative corrections. In the renormalizable field theory the latter are factorized into the 
renormalization constants of the corresponding Green functions. However, owing to the specific 
features of the renormalization in the light-like gauge, extensively reviewed in the next section, 
there is mixing of correlation functions due to the renormalization of field operators. The latter 
fact is closely related to the matrix nature of renormalization constants of the elementary Green 
functions in the axial gauge. For example, after renormalization of the fermionic propagator the 
matrix structure of the bare cut vertex could be changed, in general, since the renormalization 
matrix acts on the spinor indices of the vertex (see eq. (|37|)). 

In the Leading Logarithmic Approximation (LLA) there is a strong ordering ]TB| , |T7| ] of 
transverse particle momenta as well as their minus components, so that only the particles 
entering the divergent virtual block S, T, II or the particles adjacent to the cut vertex can 
achieve the maximum values of \k±\ and a: 

\k n ±\ < ••• < |A;2i.| < 1*1X1 < A 2 , 

a n <C ... <ti a 2 <C ai, (18) 
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while the plus components of the parton momenta are of the same order of magnitude 



~ z 2 



(19) 



for n-rank ladder type diagram. 

The radiative corrections to the bare cut vertex can be calculated using the conventional 
Feynman rules with the following modifications: 

• All propagators and vertices on the RHS of the cut are hermitian conjugated to that on 



• Every time crossing the cut have the propagator 1 / (k 2 — m 2 + zO) replaced by —2iri5(k 2 — 



• For each propagator crossing the cut there is a ^-function specifying that the energy flow 
from the LHS to the RHS is positive. (In the infinite momentum frame this is the plus 
component of the four- momentum.) 

These statements complete the rules to handle the cut vertices. 

4 Abelian evolution 

In this section we show our machinery on a simple example of abelian evolution and generalize 
it afterwards to the Yang-Mills theory. We start with overfull set of the cut vertices defined 
by eqs. (|7|)- ([T2|) and disregard for a moment the relation between them. Then eq. (|13]) verifies 
that the evolution equations thus obtained are indeed correct. We have to note that since 
the observed particle h is always in the final state some cuts of Feynman diagrams are not 
allowed and, therefore, we could not obtain the evolution kernel for the cut vertex taking the 
discontinuity of uncut graph as we are restricted over the limited set of the cuts. 

4.1 Renormalization in the light-cone gauge. 

The peculiar feature of the light-cone gauge manifest itself in the existence of additional UV 
divergences of the Feynman graphs which are absent in the usual isotropic gauges. They appear 
due to particular form of the gluon propagator 



the LHS. 




d, 




(20) 
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that possesses an additional power of the transverse momentum k± in the numerator. For 
our practical aims we limit ourselves with the calculation of the one-loop expressions for the 
propagators and vertex functions. This is sufficient for reconstruction of the equations in the 
LLA using the renormalization group invariance. 

The unrenormalized fermion Green function is given by the expression 

G- 1 (k)=f6-m -'E(k), (21) 

where E(A;) is a self-energy operator. Calculating the latter to the one- loop accuracy with a 



principal value (PV) prescription for the auxiliary pole in the gluon propagator in light-like 
gauge we get the following result 

G(k) = (1 - EOL/^AO-jr^— EM*), (22) 

f — m 



where 



777 I 

Ui(k) = 1 - — £ 2 (fc) 7+ - 7r(S 2 (A:) - XihJ, (23) 

771 I 

U 2 (k) = 1 + T-S 2 (fc) 7+ + — (E 2 (fc) - E x ) ^ 7+ , (24) 



and 



S 1 = £lnA 2 , S 2 (A;) = |-lnA 2 /^ 7 -^ 7y e? 1 (.',.'-,). (25) 
Here m is a renormalized fermion mass related to the bare quantity by the well known relation 

m = m(l-3Ei). (26) 

The functions 0™j 2 j and some useful relations between them are written down explicitly in 
the appendix. The renormalization constants are not numbers any more but matrices acting 
on the spinor indices of fermion field operators. 

An abelian Ward identity leads to the equality of the renormalization constants of the gauge 
boson wave-function and a charge Z 3 = Z g , so, that the corresponding logarithmic dependence 
on the UV cut-off cancels in their sum in the evolution equation and, therefore, we can neglect 
the fermion loop insertions into the boson line (in the QCD case this is no longer true). 

One can easily calculate the vertex function to the same accuracy. The result is 

T p (h, k 2 ) = (1 + S 1 )[/r 1 (^i)^E/ 2 (A; 2 ), (27) 

where 

Q P = 1 P - (fa - m)Q p (k 1 , k 2 ) 1+ - j+Qpifa, k 2 ){^ 2 -m) (28) 
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and 

Q P (h, h) = £ 3 (fci)7-7+7p + £3(^)7^7+7-, (29) 

here 

E 3 (fci) = In A 2 / d /^'/^ Glniz', z' - Zl , z' - z 2 ). (30) 

Apart from the graphs we are accounted for there exists an additional UV divergence of the 
virtual Compton scattering amplitude, however, we do not need its explicit expression for our 
practical purposes. This completes the consideration of virtual corrections which cause the 
logarithmic dependence on the UV momentum cut-off of the quantities in question. 



4.2 Sample calculation of the evolution kernels 



As we have noted above the UV divergences also occur in the transverse-momentum integrals of 
partons interacting with a bare cut vertex. To extract this dependence properly it is sufficient to 
separate the perturbative loop from correlation function in questions. To this end the latter can 
be represented in the form of momentum integral in which the integration over the fractional 
energies of the particles attached to the vertex is removed 



AC) 
M(0 



d A k 
W) 4 



5(( - z) 



-7717+ 



F(k), 



where 



F(k) = J d 4 xe lkx (0\^{x)\h,X)(h,X\^{0)\0). 



(31) 



(32) 



In the same way we can easily write corresponding expressions for the three-particle corre- 
lation functions. 

Let us consider, for definiteness, the fracture function X. Simple calculation of the one-loop 
diagram for the 2^2 transition in the LLA gives 



AC) 



A 2 



d ' k F{k ) f pl S (C - z")6{z" - z) 6({k " ~ k)2) 



(2tt) 4 v 'J (2tt) 3 
x{-d^(k" - k) ltl (f + m)I{f + m) lv ) 



k" A 



5(z" - Q9(z" - z ) , l2in , , „ f , „ 
2tt 2 7 (2tt) 4 ^ v "' J (z"-z) 



4 k r 



X 



a f d^k 

( 
1 

[2a" z" + kf ] 2 
dz 



d 2 k'[ / da" 5 a" + 



k 



112 



2(z"-z) / 



[2a" z" + kf ] - 2m 7+ 



a „ | [2a"z" + kf 



l(z)-M(z) 1 + 2 



{C-z). 



(33) 
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As long as logarithmic contribution appears when |fc±|/|fc"| <C 1 and a jo" <C 1 we expand 
the integrand in powers of these ratios keeping the terms that do produce the logarithmic 
divergence. Similarly one can evaluate the transition amplitudes of X to the three-particle 
correlation functions Zf 



X(C)A 2 



a 
a 



lnA 2 y dzdz'e(C-z)Z l (z',z) 



2 1 

+ 



2tt 



— — In A / dzdz'6((- z)Z 2 (z,z') 



(C-z) (z-z'Y 

(C-z) 
(C-z')(z-z'Y 



(34) 
(35) 



Due to the non-quasi-partonic |14| form of the vertex / there exists an additional contribution 
to the evolution equations coming from the contact terms resulting from the cancellation of the 
propagator adjacent to the quark-gluon and bare cut vertices. As an output the vertex acquires 
the three-particle piece 

d A k d A k' 



AO 



A 2 



(2tt) 4 (2tt) 



rZ lp (k', k)igT p (k - k', k)tG(k)5(( -z) + (c.c) 



-|-lnA 2 |^%(z',C)/ C(C * Z " ) &Uz",z"-(^"-C + z'). (36) 



As can be seen eqs. ( p3|) and (|34|) possess the IR divergences at z — £. They disappear after we 
account for the virtual radiative corrections (renormalization of the field operators) discussed 
in the previous subsection. The net result looks like 



t r = (i-e 1 )u 1 tu 2 ~ 1 , r= (/, ~ ? 



,_L_ 



: m7+, g7*7+ j • (37) 
Assembling all these contributions we come to the evolution equation for X given below by 



eq. 



4.3 Evolution equations. 

Now following the procedure just described it is not difficult to construct the closed set of the 
evolution equations 



M(C) = [-6((-z)P mm (±)m(z) 

Z7T J Z \Z j 



(3f 



dz' 



z) 



C A z(c 

-zV Zl{ ^ Z) - (C,-z<)(z-z<) 



M(z) 
2 2 (z,z'] 



(39) 
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a 

2n I 



+ 0(0 

dz 



1 



(C - CO 

c 



1(0 -M(() 



( ^_ c -X(C-C , )-^(C-C / ) 



+ 
+ 



z 
dz 



p 



'c cf 



z(C - z) s 



cc^-c+co 



eSn (C, C - C, C - C + z>) {C C c , + z,) Zi (z', 



C(*'-C + C0' 



C(c-zO(C'-^)' 



,(40) 



^(CCO 



(C'-C) 



^^(C-CO ^777^X(C0 - 7 ^(0 



cc 



c 



dZ 9(C-z)^ 7 ^-Z 1 (z-Cz) 



+ 



2 



o(C - CO (C ~ CKC ~ C + W , CO 



Qz 



A) , I I . I , C 



e^cOC-*' 



(C - *0 



-[z 2 (c,zO-^(C,CO] 



+ -2 2 (c,C0^ 



(41) 



where the dot denotes the derivative with respect to the UV cutoff ' = A 2 d/dA 2 and splitting 
functions are given by the following equations 

2 



Pmm ( z ) 
Pti(z) 
Pjm(z) 

Piz 1 (z,y) 

Pz 1 z 1 (z,y) 



z(l - z) 



+ - + 1, 

z 



-1 + 2^-1), 



z(l - z) 
2 

z(l ~ z) 
2 

z(l - z) 



2 

+ - + 

z 

+ - + T^— -S(z-l)-\n(l-y), 
z 1 - yz y 



2 

+ - + 



y 



z 1 — yz 



+ S(z-1) 



ln(l - y) 



(42) 

(43) 
(44) 

(45) 

(46) 



Now it is an easy task to verify the fulfilment of the equation of motion ( pT3|) for the correla- 
tion functions as a consistency check of our calculations. By exploiting this relation we exclude 
X from the above set of functions and reduce the system to the basis of independent gauge 
invariant quantities {Ai, T^j}- 

An important note is in order now. As distinguished from the DIS case the above eq. 
(EH) has the logarithmic dependence on the ratio of the parton momentum fractions. The 
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consequence of its presence is obvious. Taking into account the restrictions imposed by eq. 
(|17|) we can define the moments of the correlation functions in the following way 

M n = l° C f;M(0, (47) 

tc = /°°^ J o ( dcc m nc,o. (48) 

We find for two-particle cut vertex 

a 

M n = — [-i>{n + 1) - V(n - 1) - 2 lE \ M n . (49) 

where je is a Euler-Masceroni constant and ip{n) = T'(n)/T(n). And we see the universality of 
the evolution kernels for the time- and space-like quasi-partonic two-particles cut vertices, i. e. 
the Gribov-Lipatov reciprocity relation is fulfilled [|17[]. However, it is impossible to write down 
the finite system of equations for any moment of the three-parton correlation functions as the 
logarithm of the ratio of the parton momentum fractions in the evolution kernels leads to the 
infinite series of the moments as distinguished from the deep inelastic scattering where the rang 
of the anomalous dimension matrix was finite and increases as the moment of the correlation 
function increases. Therefore, it is not possible to solve the system of equations successively 
in terms of moments as well as we do not succeed in solving it analytically in a general form. 
However, in the next section when dealing with the QCD evolution we will find that the system 
can be reduced to the single equation in the limit of infinite number of colours and there is 
believe that its solution can be found analytically. 



5 Non-abelian evolution. 

For the non-abelian gauge theory the equality of the renormalization constants Z g = Z% no 
longer holds, so, we should account for the renormalization of the gluon wave-function as well 
as for the renormalization of charge explicitly. For this purposes to complete the renormaliza- 
tion program outlined in the preceding section we evaluate the gluon propagator to the same 
accuracy. The result can be written in the compact form 

D» v (k) = (l + n*(*)) U, P (k)^^U au (k), (50) 

where 

U, v {k) = g, u - ln^(fc) ^ + ^ (51) 

and 
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n- w = IhA-a/d, ^" - « + ^ c)2 + 2C41 e; lfe , - C) (52) 



are the transverse and longitudinal pieces of polarization operator. While the renormalized 
charge is given by the well known "asymptotic freedom" formula 

'N f 11 



go = g 



1 + — In A 2 

47T 



C A 

3 6 



(53) 



The second addendum arises from the diagrams with triple-boson interaction vertex. Gath- 
ering these contributions together with equations obtained in the previous section (with colour 
group factors accounted for properly), we come to the final result 



Zi ( C, 



a 



2tt\ 



C f 



0(C) ^^Mic - C) - e° n (C, C - C)^-M(C) 



+ i d io(c- Z ) 



p ZlZl ^ ^(Cz)- (c t V 



2 J CC'(z -c + o 

z z{z + CV 



(z-C + C) 



cc 



\Z 1 (z-C + C,z) 



+ J dz' 



+ c F 



-c F (q° u (C, C - O^AziO*', + eiOj^Mz', c - (') 



c, 



eSn(C, C - C, C - C + ^ C) 



C(^-C + C0 



+ ^ (eg^c, c - c, c - a ~ v ;,; / s y ztf, o 



C'(C'-^ 



^ (C + 



+ 0(0 



(('-z')(C-z' 

i (C-C' + ^)' 



c* 1 



z 1 (z'-Cx-(') 



- 20? 1 (c / ,C , -^)^7^7y^i(^O-^(C , ,C)]) }, 



(54) 



z 2 ( C, CO 



27 



-c^(c-C0^(C) 



+ /^(c-,)[p M g,g^c,.)-(c,-%) 



+ ^(-2z^l dvZ 2 (z,('-v((-z))l 



+ J dz' 



d r 1 
d^Jo 

C F 9(C-C)^Zi(z',C) 



e 



>-c+co c 



— 2 — 2 2 (^,2-C + C) 
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- e° n (C, C - / )^7y [^(C, A - MC, CO] 

- &°n(C - C, C - A - MC, (')}) 



2 1(3 U (C-*') WC'-*') 1 c 



1 * + £W.c' 



where 



(55) 



Pz 1 z 1 (z, y) = C F 



+ C 



r 


2 




|z(l-*)J 



+ - + *(* - 1) 



2 



1 - yz' 



22^2 



Ca 
2 





4 


4 


{- 




+ - 






+ z 



ln(l - y) 



l + z 



1-yz 



+ 8{z-l) 



~CV-^(lny + ln(l-y)) 



(56) 



These equations should be supplemented by the equation for the mass cut vertex M. which 



differs from its abelian analogue (38) only by the group factor Cp. 

One can easily observe the significant simplification of the above evolution equations in the 
limit N c —>■ oo since Z 2 decouples from the evolution equation for Z\. Therefore, discarding 
the quark mass cut vertex we obtain homogeneous equation which governs the Q 2 -dependence 
of the three-parton correlation function Z\. Similar phenomenon has been found in the evo- 
lution equations for chiral-even and -odd distribution functions in DIS ]10[ pl| where it has 
been observed that for multicolour QCD the momentum fraction carried by gluon in matrix 
element of quark-gluon operator varies only among the quarks ones and does not exceed the 
latter. Owing to this feature the solution of approximate equations becomes possible. In the 
forthcoming paper we will address to this issue in greater details. 



6 Conclusion. 

We have developed the effective technique for construction of the evolution equations for the 
twist-3 cut vertices. The physically transparent picture which appears in the light-cone gauge 
(an essential ingredient of our method) makes the calculations simple. Using this technique we 
have constructed the basis of chiral-odd non-polarized cut vertices closed under renormalization 
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group evolution. The identity provided by the equation of motion for field operators makes the 
construction of the basis of independent operators trivial. 

The most striking difference from the DIS evolution is the appearance of the logarithmic 
dependence on the parton momentum fractions in the evolution kernels which makes the suc- 
cessive solution of equations in terms of moments impossible. The nonlocality of the cut vertex 
in coordinate space is essential since even if we start from the local cut vertex it will smeared 
along the light cone upon the renormalization. 

However, an important point is that there is decoupling of three parton correlation function 
Z\ from Z 2 . The situation has the closer similarity with DIS where in the limit N c — > oo 
there was a very important simplification as the evolution kernels have been vanishing for 
contributions with interchanged order of partons on the light cone, i.e. the gluon momentum 
fraction ranges between those of the quarks only. This property allowed the author of refs. 
|i~0| , |l~9f to find the solution of simplified equations exactly on the nonlocal form. In our case 
the decoupling of Z\ may have the same consequences. This work is in progress now and the 
results will be published elsewhere. 

We are grateful to L.N. Lipatov for helpful discussions and encouraging. One of the authors 
(A.B.) would like to thank CERN Theory Division for hospitality extended to him during his 
visit where this work was finished. He benefited from useful discussion with V.M. Braun, D. 
Miiller and M.A. Shifman. This work was supported by Russian Foundation for Fundamental 
Research, grant N 96-02-17631. 



A Appendix 



In this appendix we present the necessary properties of the O-functions 




(57) 



For our discussion it is enough to have an explicit form of the function 



©11(21,22) 




(58) 



z\ - z 2 



since the others can be expressed in its term via relations 




(59) 



zi - z 2 




(60) 



(61) 
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In the main text we have used the identity 

c 



PV / dz 



'(C-z) 



e u n (z,z-Q+& n (C,C-z) =0. (62) 
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